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SUMMARY 


A procedure that allows preliminary assessment of the effects of component size 
changes on the characteristics and performance capability of aircraft has been 
developed into a computer program. Applications can be made to subsonic, super- 
sonic, and hypersonic aircraft using JP, liquid-hydrogen, or liquid-methane fuel. 
Jet engines as well as propeller systems can be treated. Basically, the program 
requires an input baseline aircraft on which component sizing changes are made; new 
aerodynamic, propulsion, and weight characteristics are then determined and a 
mission profile, including a reserve fuel segment, is computed to find the new 
aircraft's range capability. One of the output options provided by the program is 

the basis for preparing aircraft sizing "thumbprints " which pinpoint the wing and 
engine size for best performance. 


INTRODUCTION 

A computer program for aircraft sizing is a useful tool for preliminary aircraft 
system studies. The method must account for the interactions between aircraft 
geometry, aerodynamics, propulsion, and weights while providing rapid, low-cost 
evaluations of design option:,. Important advantages offered by this capability 
include the realistic assessment of technological advances and the identification of 
promising areas that justify further stjdies in greater depth. 

A program that satisfies these requirements has been developed for use in the 
Aeronautical Systems Office at the NASA Langley Research Center. An c-arlier version 
was reported in reference 1. Briefly, the method performs sizing operations on a 
baseline aircraft and then a mission analysis to determine the effects of sizing on 
range. In the few instances where comparisons are possible, results of the method 
are in close agreement with similar ones in industry (refs. 2 and 3). 

This paper rnntains * description of the processes used in the oro>jr?m, some 
applications of the program to subsonic, supersonic, and hypersonic aircraft, and 
definitions of Input and output. For identification purposes, the program is called 
ASP, (A)ircraft (S)izing (P)erformance program. 
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Subscripts 


AP 

approach 

BL 

baseline aircraft or bleed drag 

F 

fuel 

LO 

lift off 

MAX 

maximum 

SL 

sea level 


A dot over a symbol denotes its time derivative. 


PRELIMINARY CONSIDERATIONS 

The conputer program described herein cannot be used to design an aircraft 
concept. Instead, the program performs operations on a baseline aircraft which has 
already been designed in sufficient detail to provide the required inputs. Within 
the program, this baseline aircraft is sized according to optional input values such 
as wing loading, engine thrust, number of passengers, and gross weight. After 
sizing, the aerodynamic, propulsion, and weight characteristics of the aircraft are 
determined; then a mission profile with reserve fuel segments is computed to find 
the range of the aircraft. 

Some of the output provided by the program include: 

1. Range for given gross weight and number of passengers 

Z. Gross weight for given range and number of passengers 

1. Number of passengers for given range and gross weight 

4. Payload - range curve 

5. Passenger and fuel offloading 

6. Effects of technology advances in aerodynamics, propulsion, and weights 

7. Effects of in-flight refueling 

8. Mission radius to payload drop 

9. Information for preparing an aircraft sizing "thumbprint.” (A "thumbprint" 
is a map of constant value contours, like range or gross weight, plotted as 
a function of engine and wing size.) 

The program can he applied to subsonic, supersonic, or hypersonic aircraft using 
JP, liquid-methane, or liquid-hydrogen fuel. Jet or propeller propulsion systems 
can be considered. 
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Five baseline aircraft have been chosen to Illustrate the capabilities of the 
program. Computer drawings of these aircraft are shown in figure 1. All aircraft 
are preliminary concepts designed during in-house activities. 


METHODS 

Component Sizing 

Fuselage .- To simplify drag calculations .nd packaging, the fuselage is 
approximated by an equivalent body consisting of a constant area midsection with end 
caps having a Sears-Haack profile (ref. 4). The constant area midsection contains 
the passengers, door-gal ley -lavatory sections, cargo, and, if required, cryogenic 
fuel tanks. Variables that affect the passenger compartment are the number of seats 
abreast, seat pitch and width, aisle width, and the length and number of door- 
galley-lavatory sections. Although these quantities may be input, their values may 
not all be preserved because final values are determined from the baseline aircraft 
using input values for passenger compartment length and abreast seating arrangement. 

The arrangements available for fus' ' :ge packaging are shown in figure 2. The 
first arrangement (fig. 2a) applies to aircraft witn one-passenger level and conven- 
tional fuel. In addition to the passengers, the arrangement contains the cargo 
normally stored benea h the floor; however, if not enough space is available, the 
excess cargo is stored in a section behind the passenger compartment. The cargo is 
assumed to be stored in cargo containers which occupy 60 percent of the available 
volume. 

If the aircraft uses cryogenic fuel, the second arrangement (fig. 2b) places all 
the cargo behind the passenger compartment and fuel tanks are located beneath the 
floor. If more fuel volume is required, it is located in full-depth sections of the 
same length at both ends of the midsection. 

The third arrangement (fig. 2c) has two identical passenger levels with aft 
cargo storage. If cryogenic fuel is required, it is stored in duplicate tanks fore 
and aft of the passenger compartment. 

In application, the midsection length and diameter are calculated first for the 
requirements of the baseline aircraft. Then, from the length of the baseline 
aircraft, the length of the fuselage end caps are found. During sizing, a new 
midsection of different length and diameter is determined to satisfy the require- 
ments for passengers, cargo, abreast seating, and fuel volume (if the fuel is 
cryogenic). The new fuselage length is then obtained by adding the end caps to this 
midsection. 

Wing and Tail Surfaces .- From the Inputs for the baseline wing and tail 
surfaces, definitive geometric properties are determined and then normalized by a 
representative dimension (such as wing root chord or geometric chord). Baseline 
tall volume coefficients are also found for each surface. During sizing, the non- 
dimensional properties are used to shape the components In a geometrically similar 
manner. Wing size is specified by wing loading or reference area, while the tail 
surfaces are sized to keep their volume coefficients equal to baseline values. The 
locations of the components are established by keeping mean aerodynamic chord 
positions (25-percent subsonic, 50-percent supersonic) at the same percent body- 
length station as the baseline aircraft. If overlap occurs between the wing and 
tail surfaces, the program will output a warning and continue operation. 


ORIGINAL 
OF POOR 

Fuel Volume . - Conventional fuel Is assumed to be stored In the wings and in the 
fuselage where wing carry-through structure Is located. The maximum wing-fuel 
volume, V^, is given by 


pAGEtf 

QUALITY 


V V W,EXP = K A + % V w',tXP 


where V w ^ Is the total volume of the exposed wing, a constant determined 

from the baseline value of maximum wing-fuel volume, and = 7,16 X 10~®, which 

was derived from studies of wing-fuel capacity performed by the Kentron 
International Inc., Kentron Technical Center, Hampton, Virginia. The maximum 
fuselage-fuel volume Is scaled from the oasellne value and it is assumed to vary 
directly with root chord, maximum thickness of the wing, and the maximum width of 
the fuselage. 

Cryogenic fuel is contained only In the fuselage at locations covered previously 
under fuselage sizing. Fuel tanks are simple constant area tanks having a length 
that provides the required volume. The tank cross-sectional shape Is the same as 
the fuselage after accounting for the wall thickness for cryogenic tankage. 

Propulsion System. - The length, L,, and diameter, 0 N , of circular nacelles are 
sized by the Installed, sea level, static thrust, T $L , In the following manner: 

N l 

L N = L N,bJ T SL /T SI., 8|J 
°N = J N,8L- T Siy T S!.,BL^ 


a 


where the subscript 31 refers to the baseline aircraft and T $L varies with Input 

thrust-to-welght ratio (or power-to-welght ratio for propel ler-dri ven aircraft). 
The engine sizino exponent N|_, which varies with engine concept, is an input. 

For two-dimensional inlets, the wetted area is assumed to scale directly with 
thrust. Both the length and width of these Inlets, therefore, vary directly with 

^SL^SL.BL^ 5 * 

In a similar manner, the propeller diameter, Dp, is sized by 


D P s D P,8L‘ T SL /T SI.,BL^ 


.5 


f 



and the RPM of the propeller, RPM p , and engine, RPM £ , by 
RPHp - «PMp, BL (T SL>BL /T SL )- 5 

RWI E = rpm EiBL ( t sl, B u /T Sl )' 5 

Nacelle positions are established by keeping their lengthwise locations a constant 
percentage of the local wing chord and their spanwlse pod-to-pod and pod-to-fuselage 
clearances a constant percentage of wing span. 

An Illustration of the effects of these sizing procedures on the supersonic 
transport Is shown In figure 3. 


AERODYNAMICS 


The baseline aerodynamic Input Includes lift and drag polars for both takeoff 
and higher speed conditions and zero-lift drag Items. All, except the takeoff data, 
are corrected for sizing effects. Corrections are obtained by applying the methods 
described below to both the baseline and the sized aircraft, then adding the 
difference in the results to the baseline input value. 

Llft-Oraq Polars. - For subsonic aircraft, the shape of the lift-drag polars are 
functions ot wing aspect ratio, thickness, and sweep, and also Mach number when 
compressibility effects become important. To account for these variables, the 
program uses a method based on the drag build-up system given in reference 5. Some 
results that the method provides for sizing variations to the subsonic transport are 
shown In figure 4. 


For higher speed aircraft, the program treats only wings that are geometrically 
similar to the baseline (equal aspect ratio, sweep and thickness ratio). It is 
therefore assumed that the shape of the lift-drag polars for these aircraft is not 
affected by changes In component size. 

The lift-drag polars at all speeds will shift along the drag axis because of 
drag increments that result from sizing effects on zero-lift drag. 


Zero-Lift Vag. - The Items that contribute to 


Cj^ are the friction drag, 


roughness drag, and wave drag. To these, engine bleed and air-conditioning drag 
must be added If they are not Included In the baseline engine data. All drag items 
are corrected for sizing effects. The drag coefficients for roughness, C D RUF , air- 
conditioning, C n ac , and engine bleed, B , vary by 


C 0,RUF = ^ C D ruf ^ bl ^ S WET /S REF^^ S WET /S REF^BL 
C n,AC s ( c o,acV ( t sl /s ref^ t sl /s ref)bl 
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C D t B = ^D.Bfo. ^ T SL /S REF^^ T SL /S REF^BL 

where Syr-r is the aircraft wetted area, NBAS, the number of passengers, and S RE p, 
the aircraft reference area. 

Wave Drag .- To avoid the long computing times required to determine the wave 
drag by existing computer codes, a quick, simple, but approximate approach was 
adopted. The method treats the aircraft as a collection of isolated components, all 
t exposed to free-stream conditions. Since the wing and tail surface vary geometri- 

cally along with the aircraft reference area, the wave-drag coefficients for these 
surfaces are assumed to remain constant as the size of the component changes. 

Sizing, however, does affect the fuselage wave-drag coefficient. As mentioned 
» previously, the fuselage is composed of a cylindrical midsection with identical end 

caps. Only the end caps are assumed to contribute to the wave drag and their shape 
was chosen to be a Sears-Haack profile which allows the fuselage wave-drag coeffi- 
cient to be obtained (see ref. 6) by 

C DW,F = 24V ( S MAX /S REf) /L 

where V, l, and S MAX are the total volume, length, and maximum cross-sectional 
area of the fore- awT^ft-end cap combination. 

The contribution of the engine nacelles to the wave drag is considered to be 
small enough to ignore the effects of sizing, 

A comparison of the results of this simple approach with those from the wave- 
drag program of reference 7 for the supersonic transport aircraft with JP fuel is 
shown in figure 5. Although some differences do occur, the method provides an 
improvement over the assumption that wave drag is unaffected by sizing. 

Friction-Drag.- For calculating skin-friction drag, the fuselage, tail surfaces, 
and engine nacelles are treated as flat plates with corresponding wetted areas and 
lengths. The wing, however, is treated as a collection of strips with the strip 
mean geometric chord as the characteristic length. All surfaces are at zero angle 
of attack. Skin-friction drag corrections are made at all Mach number-altitude 
combinations required by the mission flight profile. 

Two different procedures are used to calculate friction coefficients. One, 
which applies to fully turbulent flow, uses the reference temperature method 

* (ref. 8) along with the constants for turbulent flow given in reference 9, and the 
Prandtl-Schlichting formula for incompressible average skin friction (ref. 10). An 
input value for emissivlty is required for the solution of wall temperature. 

* Tho other procedure applies to mixed laminar-turbulent flow and is based on 
assumptions of an adiabatic wall temperature, a recovery factor of .88, and a dis- 
continuous change from laminar to turbulent flow at the specified location. Since 
these assumptions are more appropriate to wind tunnel, rather than flight 
conditions, the results are not used directly. They are used, Instead, to Increment 
the fully turbulent flow results obtained as described above. 
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An illustration of the effects of sizing on the zero-lift wave drag and the 
maximum lift-drag ratio for the supersonic JP-fueled transport under turbulent flow 
conditions is shown in figure 6. 


WEIGHTS 

The weight input for the baseline aircraft include component structural weights, 
propulsion system weights, system and equipment weights, and crew, passenger, and 
cargo weights. Empirical weight relationships are used to find the effects of 
sizing on each item. These emperical weight equations and the method in which they 
are used are described in Appendix A. 

An example of the weight results obtained from sizing the supersonic JP-fueled 
transport is shown in figure 7. 

MISSION ANALYSIS 

After the aircraft is sized and its aerodynamic, propulsion, and weight 
characteristics are known, a mission analysis is made to determine the range 
capability. The mission profile includes taxi, takeoff, cl imb-to-cruise Mach 
number, climb-to-cruise altitude, cruise, descent, approach, and reserve segments. 
For supersonic cruise missions, a subsonic cruise leg may also be specified for the 
outbound and/or inbound flight segment. 

Taxi and Takeoff.- During taxi, the engines a-e throttled to provide the fuel 
flow rate and thrust required to overcome rolling friction. The fuel for taxi is 
found from this fuel flow rate and the input taxi time. 

The Mach number-altitude profile during takeoff is not input; it results from 
Influence of the aircraft's aerodynamic and maximum power characteristics on the 
equations of motion. 

^ - (T cos a - 9) /VI - sin 9 

9 = (T sin a + L)/W - cos 0 

= V sin 9 

1$ ■ V cos 9 

v? = - vJp 

During the ground run, these equations can be solved easily because 9 and 
are zero. But after "lift-off," when these simp! i fications no longer apply, they 
are solved by more time-consuming numerical integration using the fourth order 
Runge-Kutta formula. 

The takeoff segment Includes the ground run and climb over the obstacle to an 
altitude of 400 feet. Analyses Include '“all-englnes-go" ar.d "one-engine-out" 
situations but no analysis of the balanced takeoff field length (ref. 11) is made. 


As shown in figure 8, the the ground run begins at zero velocity, zero flap 
deflection, and ground angle of attack using takeoff aerodynamics that include 
landing gear drag increments and ground effects. At velocity ^rotate* the flaps 

are deflected and the aircraft begins its rotation to the takeoff angle of attack. 
After rotation is complete, the angle of attack and flap deflection are held 
constant until "lift-off" occurs. After "lift-off," the aircraft climbs with angle 
of attack and flap deflection unchanged; however, the aerodynamics are corrected for 
landing-gear drag retraction over an input time interval and for the diminishing 
influence of ground effect over an Input altitude interval. The climb continues at 
constant angle of attack until the input climb gradient is reached; then, angle-of- 
attack changes are made to maintain this gradient. The end of the takeoff segment 
occurs when the aircraft reaches an altitude of 400 feet. 

When the aircraft reaches the height of the obstacle (typically 35 feet) 
115 percent of the horizontal distance covered is the takeoff field length when all 
engines are operating; but, for the "one-engine-out" condition the takeoff field 

length is the actual distance. 

The minimum takeoff field length, as indicated in reference lx, occurs at an 
optimum value of the velocity for start of rotation, V ro t a te* If more than a few 

aircraft are to be analyzed, however, the computing time required to find the 

variation In takeoff field length with V ro tate becomes prohibitive. A shorter 

approach was therefore adopted in which an Iterative solution Is used during the 

ground run to find the value of V rotate that allows rotation to be completed and 

the takeoff angle of attack to be attained just at "lift-off." The takeoff field 
length that results from this value of V rotate is a close approximation to the 
minimum value as shown in figure 9. 


Cl i mb . - After the takeoff phase is complete, the aircraft is at an altitude of 
400 feet. Its Mach number and altitude are then phased into the input H-M profile 

and the climb segment begins. This and all remaining segments of the flight profile 
use simplified differential equations of motion (ref. 12). Theie simplifications 
are based on the following assumptions: the time derivative of the flight-path 

angle is zero, the angle of attack and flight-path angle are small, and the 
aircraft'? altitude is negligibly small compared to the earth's radius. For climb, 
over Interval points 1 and 2 along the specified flight path, the weight change, 
when velocity is increasing, is given by 


In (W 2 /Wj) 


-5TT \l ? - V t + g (DH/DV) In (Vg/Vj) 

__ (V- 3/T) ~~~ 


and, where velocity is constant, by 


In (W 2 /Wj) = 


■SKI. « 2 - 

V 1 - 0/T 


OH/DV is the climb-path slope and the barred quantities represent effective values 
which are taken to be simple averages of Interval point values. 
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Time and range changes over the interval are: 




(», * v zK 

*2 -*i- -t — - t 1 ;- *i ) 

where iJp is the average fuel flow over the interval. 

The aircraft follows the input H-M profile up to the cruise Mach number unless, 
for a supersonic cruise mission, an outboard subsonic cruise leg is specified. If 
this occurs, the climb stops at the input subsonic cruise Mach number and a search 
(at constant M) is made for the altitude that provides the best Breguet factor with 
engines throttled. The subsonic cruise leg is then flown at constant lift 
coefficient for the specified range. Next, the aircraft accelerates at constant 
altitude until it intersects the input H-M profile, which it then follows to the 
supersonic cruise Mach number. During the climb, the acceleration can be controlled 
by input to occur at maximum power or at partial power settings that use least fuel. 

After the cruise Mach number is reached, a climb at constant Mach number is made 
to the cruise altitude. Unless specified otherwise, cruise will start at the 
altitude which provides the best Breguet factor; however, the altitude at the start 
of cruise can be input directly or restricted by inputs for constant angle of 
attack, constant throttle setting, or minimum rate of climb. After the cruise 
altitude is reached, the climb fuel is known. 

Reserve s.- Since the fuel available for cruise and descent is the remainder 
after taxi, takeoff, climb, and reserve fuel are allotted, the reserve fuel require- 
ments must be determined next. Reserves can include a trip fuel allowance, fuel for 
a missed approach, a flight to an alternate airport, and a hold at constant Mach 
number and altitude. A typical reserve flight profile is shown in figure 10. 

The solution for reserve fuel, which is iterative since tne weight of the air- 
craft at the end of descent is not known at this time, starts with the trip fuel 
allowance. This fuel allowance for subsonic aircraft is obtained from the fuel 
burned for an input time interval using the fuel flow rate at the start of cruise. 
For other aircraft, the allowance is on input percent of the trip fuel, Wp j, which 
at this time may be obtained from ' ’ 


W F,T “ W F * W F,R " W F,TAXI 

where Wp is the total fuel weight, Wp >TAXI the taxi fuel weight, and Wp >R the 

total reserve fuel weight (which must be assumed for the first iteration). The fuel 
for missed approach is obtained by either calculating a "wave-off" and acceleration 
to an altitude of 700 feet with landing gear retracted and "out-of -ground" aero- 
dynamics In effect, or, more simply, by the fuel burred for an input time Interval 
(typically 2 minutes) at the takeoff fuel-flow rate. For the detailed calculations, 
the weight at the start of the missed approach segment is: 
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, the trip fuel allowance. 


W * M 6 - W p + W p ^ R - W FjT p A 
where Wg is the gross weight and W R TpA 

For the alternate airport requirement, the aircraft climbs a linear H-M profile 
from the end of the missed approach to the altitude at the start of cruise, cruises 
the required distance, then descends to approach conditions. The total range for 
these segments is the required distance to the alternate airport. The analysis is 
performed for various climb and cruise throttle settings and cruise Mach number- 
altitude combinations to find conditions for least fuel. 

During the descent to the alternate airport, the hold segment is performed. The 
time for hold is input and the Mach number-altitude combination is again found for 
the least fuel requirement. The range covered during hold is not credited to the 
alternate airport distance. 

After all reserve segments are completed, the total reserve fuel is updated and 
the process repeated until successive values are within .1 percent. An example of 
the effects of sizing on the reserve fuel requirements for the superson’'. transport 
are shown in figure 11. 

Cruise .- The total fuel available for cruise and descent is now known, and 
solutions of these segments are repeated until the required and available fuel are 
within .1 percent. The process starts by assuming the fuel weight for descent 
(about .5 percent of the aircraft gross weight) then finding the fuel for cruise. 
During cruise, toe aircraft flies at the selected cruise option (like constant C L 
or others described below) until the cruise fuel is exhausted. 

For cruise, the increment in range, R, over interval points 1 and 2 is given by 


R ? - Rj * BF In [1/(1 - Wp/Wj) ] 

where Up is the fuel weight used over the interval, and IF is the numerical 
average of the Breguet factor, which changes over the interval because of the input 
cruise option selected. The computing interval during the main cruise segment is a 
fuel weight increment obtained from the weight of cruise fuel and an input number of 
intervals. For outbound or inbound subsonic cruise legs, the computing interval is 
the entire cruise leg. 

The main cruise segment is normal 1> controlled by different input options. 
These include: cruise at maximum available Breguet factor, constant lift 
coefficient (or constant angle of attack), constant altitude, and constant engine 
throttle setting. The different types of cruise segments, of course, affect the 
cruise altitude variation and cruise range of the aircraft. Some examples of this 
behavior for the supersonic transport for several wing sizes are shown in figure 12. 
For all cruise options, except maximum power, the start of cruise was selected to 
start at the altitude for maximum available Brequet factor. 

With the largest wing (W/S = 50 psf), the altitudes for maximum Breguet factor 
are above the climb ceiling. Therefore, unless constant altitude is specified, 
cruise, with the better range, will occur at maximum power. 
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With the median wing (W/S « 85 psf), the maximum Breguet factor at start of 
cruise occurs at a lower altitude (60,000 ft) and thrust is now sufficient to 
provide the altitude variation that meets this condition throughout cruise. For 
this wing and engine size combination, the constant Cj_ (or angle of attack) option 
is equivalent to the maximum Breguet factor option and any of these gives the best 
cruise range. 

With the smallest wing (W/S = 110 psf), the altitude for best Breguet factor at 
the start of cruise is below the climb path. A constant altitude cruise with 
decreasing therefore occurs until maximum Breguet conditions can be met and the 
best cruise range now occurs only for this option. 

From these results, regardless of thrust or altitude limitations, It is evident 
that the best cruise range will always be obtained by selecting the option for 
maximum available Breguet factor. 


Long-Endurance Cruise .- For this type of mission, the cruise segment must use 
minimum fuel and the program will find the variations in Mach number, altitude, and 
lift coefficient to provide this condition throughout cruise. If the optimum 
conditions prove to be unsatisfactc ry, there are inputs available to specify the 
cruise Mach number and/or altitude. An example of the effects of optimum and non- 
optimum conditions on endurance for size variations of the long-endurance aircraft 
are shown in figure 13. 


Descent. - Like the climb, the computing interval for the descent is a Mach 
number increment. It is defined by the difference between the Mach number for the 
main cruise segment and the Mach number at the end of descent and an input number of 
intervals. If velocity is decreasing, the change in time, t, over the interval is: 


t 


2 " 


W 1 V V 9 (DH/DV) 1n 'W 


g 


T - D 


and if velocity is constant 


t 


2 ’ 


s r 



Again, numerical averages of interval point values are used for the barred 
quantities and the altitude at point 2 is found at the lift coefficient for 
(L/D) max with Wo = W^ (a reasonable assumption since fuel expenditure during 
descent is very low). W^ is subsequently incremented by 


W 2 - Wj » '.J F (t 2 - tj) 

After cruise, the aircraft enters the descent phase by making the transition to 
flight at (L/D) m ^ x during the first computing Interval. The remainder of the 
descent is then performed at (L/D) MA v. If the aircraft Is supersonic and an Inbound 
subsonic cruise leg is chosen, the descent stops at the specified Mach number and a 
search is made for the altitude giving the best Breguet factor. A subsonic cruise 
at constant lift coefficient is then performed over an input range, after which 
descent to the start of the approach at M * .4 and H = 4000 feet is calculated. 
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Approach .- This portion of the descent extends from M s .4 and H = 4000 feet 
down to an aUitude of 370 feet. Using the equate . given previously for the 
descent, an iterative solution is required to find the aircraft weight and velocity 
at the end of the approach. The angle of attack and flap deflection for approach 
must be Input. The takeoff aerodynamics are used without ground effect but the drag 
increments for the extended landing gear are included. 

Overall Mission Characteristics .- Typical results from different segments of the 
mission analysis are now presented for several wing and engine sizes on the super- 
sonic transport with jP fuel. For all size changes, the aircraft gross weight was 
held constant and the cruise segment was specified to occur at maximum available 
Breguet factor. 

Figure 14 shows the flight profile that results from changes in wing size with 
T/W = .45. Controlled by input, the climb profile is restricted, but the cruise and 
descent profiles are different for all cases. 

Takeoff results, presented in figure 15, show that changes in the fuel required 
due to sizing are small; the spread of values is less than .25 percent of the air- 
craft gross weight. On the other hand, changing lift (W/S), or acceleration (T/W), 
capability has a large effect on takeoff field length. The velocity at "lift-off" 
is seen to be primarily a function of wing loading. 

Climb results that show excess thrust capability are given in figure 16. 
Decreasing with Mach number, the excess thrust reaches its minimum value at the end 
of the clifro, and this general trend is not affected by sizing. Changes in excess 
thrust are, of course, almost directly proportional to changes in engine size, while 
those resulting from wing changes are much less pronounced since they reflect only 
the different drag characteristics encountered during the climb. 

The climb time, fuel, and range are shown in figure 17. The effects of engine 
size are large and aircraft with low acceleration, naturally, take longer times, use 
more fuel, and cover more range to reach cruise. Wing size does not seem to be too 
important except for the larger wings (W/S = 70 and 50 psf) at lower engine sizes 
where time, fuel, and range increase while climbing to the higher cruise altitude. 

The fuel available for cruise has a significant effect on cruise range and time 
as shown in figure 18. Because of the additive affects of wing and engine size on 
total fuel and the fuel required for other flight phases, the aircraft with the 
smallest wing (W/S = 110 psf) and medium engine size (T/W = .35) has the maximum 
available fuel for cruise. The maximum cruise range, however, occurs at 
W/S * 90 psf because the slightly larger wing improves the L/D characteristics 
without an appreciable decrease in available fuel. The more efficient larger wings 
produce lower ranges because their high weight limits the available fuel. 

Descent results are given in figure 19. Fuel for descent is less than 1 percent 
of the aircraft gross weight for all size variations. Because the descent is 
essentialy a glide at (L/D) H ^ X , aircraft with the largest wings give the longest 
range. The effect of engine size is primarily due to different nacelle skin- 
friction drag Increments. 

Approach speeds (fig. 20) are primarily a function of wing size, although engine 
size also has an effect, but much smaller, because of its Influence on approach 
weight. 
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Total trip time and range are shown In figure 21. Because the longer ranges in j 

the climb and descent segments (figs. 17 and 19), the lowest engine size provides j 

the greatest overall range. 1 


APPLICATIONS 

To show some of the capabilities of the program, application will now be made to 
four different types of aircraft which Include: two supersonic transports, one 
using JP fuel, the other, liquid-hydrogen fuel; a subsonic propeller-driven 
transport; and a hypersonic aircraft using liquid-methane fuel. The baseline air- 
craft are shown in figure 1. During the discussion. Input variables may be referred 
to. For a definition of these variables, see Appendix B. 


Sizing Thumbprints 

During each computing cycle, selected output from the sizing and mission 

analysis are recorded on a file named TAPE14. This data can be used to construct 
the sizing "thumbprint"— a diagram that pinpoints the best wing and engine size of 
the air'~aft. The "thumbprint" contains contours of quantities that are to be 
optimized and limit contours of quantities that restrict the choice of wing and 
engine size. Examples of optimizing contours are range, gross weight, and seat 

miles/gallon, while takeoff field length, fuel limit, maximum rate of climb and 
thrust margin are typical limit contours. TAPE14 and an example of all the 

information it contains for contouring is described in Appendix C. 

A typical "thumbprint" for the supersonic transport, using JP-fuel , with a range 
of 4000 nautical miles is shown in figure 22(a). To obtain the data for this 

"thumbprint," the essential inputs are the required range (RNGDES), wing sizes 
(WOSTB), and engine sizes (TOWTB). The contours are for constant gross weight and 
they show that a minimum occurs at a wing loading, W/S, of approximately 90 psf and 
a thrust-to-weight ratio, T/W, of about .28. The limit contours, however, show that 
this combination is unrealistic because the wing is too small to hold the required 
fuel, and the takeoff field length is longer than 10,000 feet. If desiqn require- 
ments include a maximum lift-off speed of 200 K, a maximum approach speed of 160 K, 
a maximum takeoff field length of 10,000 feet and a minimum rate of climb during 
cruise of 300 ft/min., the choice of wing and engine size are those included within 
the boundaries provided by the limit contours for approach speed, takeoff field 
length, and rate of climb. Influenced by these restrictions, the conditions for 
lowest weight aircraft would change to about W/S = 82 psf and T/W = .32. 

Although a low gross weight is an important quantity that decreases the cost of 
the aircraft and its maintenance, it is not the only criteria for selecting wing and 
engine size. Quantities that relate to the cost of operating the aircraft may also 
be important. One of these, which could assume a major role if fuel prices 
multiply, is a high value of the seat-mi les/gallon-of-fuel irameter. The 
corresponding "thumbprint" for the supersonic transport with contours of this 
parameter is shown in figure 22(b). For this aircraft concept, the conditions for 
maximum passenger-miles/gallon are about the same as those for minimum gross weight 
with flight restrictions considered. 
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"Thumbprints" for the other aircraft are shown In figures 23 to 25. These 
Include the hydrogen -fueled transport, the subsonic transport, and the methane- 
fueled hypersonic aircraft. The "thumbprints" for the transports are at constant 
range, whereas the hypersonic "thumbprint" is at constant gross weight. To obtain 
this “thumbprint," inputs for design gross weight (WGDES), wing sizes (WOSTB), and 
engine sizes (TOWTB) are required. There are no fuel limit lines on the figures for 
the aircraft with cryogenic fuel because these aircraft are sized to hold only the 
required fuel. The symbols on the "thumbprints" define the gross weigh 1, and wing 
and engine sizes of the "design-point" aircraft that will be used to show some of 
the additional capabilities of the program. Although the capabilities discussed are 
divided among the different types of aircraft, applications, in most cases, can be 
made to all types. 


Supersonic Transport (JP fuel ) 

Payload - Range r raue. - With design-point Inputs for gross weight (WGDES), wing 
loading (WOSOES), thrust -weight ratio (TOWDES), and the payload-range trade 
requested (ICALPRC=1), data for the curve shown In figure 26 will be output. The 
computing process Involves a computing cycle at the design conditions, then 
additional cycles with payload off-loaded. Two possibilities can occur when the 
payload is off-loaded. The first occurs when the aircraft has ample fuel capacity, 
then as payload decreases, fuel is added to keep the gross weight constant. Range 

increases because of the Increased fuel mass fraction (fuel weight divided by gross 

weight). If, however, the aircraft cannot hold the added fuel, the gross weight Is 
decreased by an amount equal to the weight of the excess fuel. Again the range 

Increases (fuel -mass fraction has increased) but this time, not as rapidly. These 

two conditions are Illustrated by the solid and dashed curves in figure 27. 

Radius to ’ayload Prop .- A fallout from the payload-range analysis is the 
radi us-to-payl oad drop useful for bomber studies. To provide this output, the 
program executes two computing cycles--one with full payload and one without 
payload. The results are then used to find the point at which the outbound range 
with payload equals the Inbound range with paylcad out. The program does not 
account for thp possibility that the altitudes, with and without payload at the 
drop, may not be equal. Typical results for the drop radius compared to radii with 
and without payload are shown in figure 27. 

Fuel ift-LOdO - 'iy. - Distances below the design range can be flown with or without 
fuel off-loading. Both cases can be treated by the program. Reduced ranges without 
fuel off-loading are Initiated by inputs in RRTB and solved simply by shortening the 
cruise segment the required amount. Reduced ranges with fuel off-loading are 
obtained by Inputs in OLRTB and complete mission analyses are performed which solve 
for the required fuel. From the standpoint of fuel efficiency, fuel off-loading Is 
the better of the two because the dead weight of excess fuel Is not .arrled along. 
Figure 28 shows this result as well as the fact that for best fuel efficiency a 
redesign of the aircraft for the reduced range capability would be required. 

Subsonic cruise legs .- During a supersonic mission, subsonic cruise may be 
required during the outbound and/or inbound flight segments to avoid sonic boom 
problems. If these options are chosen, the range capability of the aircraft Is 
reduced because of the lower cruise efficiency at subsonic speeds. Typical range 
penalties that result from different subsonic cruise ranges are shown In figure 29. 
The best subsonic range Is, of course, obtained by all subsonic cruise because of 
the extra cruise fuel gained by eliminating the acceleration to cruise Mach number. 
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Passenger Load. - Abreast seating arrangements, to some extent, affect the f 

aircrafts range because changes In fuselage length and fineness ratio alter the 1 

drag and structural weight. This effect Is shown in figure 30 where the gross f 

weight required for a range of 4000 nautical miles with various passenger loads at | 

different seating arrangements are shown along with accompanying changes in fuselage 
length and fineness ratio. Four or five seats abreast are reasonable choices, but 
six-abreast seating is not competitive because of high fuselage drag. 


Supersonic Transport (Liquid-Hydrogen Fuel) 7 

u 

In this program, the low-density hydrogen fuel is located only in the 

fuselage. Since the fuel weight is affected by fuselage size and weight, which, in 
turn, changes with fuel weight, aircraft sizing is an iterative process which 
continues until the contained fuel equals the required fuel. Because of the large * 

fuel volume, the maximum equivalent cross-sectional area and length of the fuselage 
is not only dependent on passenger load and seating arrangement, as is the case for 
noncryogenic fuel, but also on the volume of the fuel tanks. In order to make the 1 

fuselage sizing a tractable problem, the process may be controlled, through inputs, 
to occur at constant fuselage length (KBL), constant maximum equivalent radius 
(KBR), or constant fuselage fineness ratio (KBFN). The effects of these choices on 
the range capability of several size variations of the aircraft are shown in 
figure 31. 

An alternative solution to passenger packaging for these large volunie aircraft 
is the use of two-passenger levels (LEVELS - 2). With this choice, the fuselage 
maximum equivalent cross-sectional area is again dependent on the passenger-abreast 
seating arrangement and the fuselage length and fineness ratio cannot be controlled 
by input. For comparison, the range for aircraft with two-passenger levels are also 
shown in figure 31. For these aircraft, two-level arrangements may be superior and 
this inference agrees with industry conclusions in reference 13. 


Subsonic Propeller Driven Transport 

The performance characteristics of this type of aircraft depend on the design 
conditions (altitude, Mach number, and tip speed) of the propeller. These 
conditions, which cannot be controlled in the program, must be used to develop the 
input characteristics of the baseline engine discussed in Appendix B. The effects 
of these propeller design conditions on the baseline aircraft and some of its size 
variations, as provided by this program, are discussed in reference 14. 

Technology Improvements .- The effects of technology advances are found by simple 
percentage changes in the required areas with interactions in other areas falling 
out as a result of the sizing and mission analysis procedures. Areas available for 
improvement (or decav) are given in Appendix B. To illustrate some results, the 
sensitivity of gross weight and seat miles/gallon for the subsonic transport to 
chafes in engine weight, propeller weight, specific fuel consumption, drag, and 
structural weight are shown in figure 32. 

Laminar-lurbulent .low .- Another technical advance that shows considerable 
promise. If successful, Tor improving the energy efficiency of the aircraft is 
laminar flow control which preserves laminar flow over large areas of the aircraft 
components. Potential gains, predicted by the program, are shown in figure 33 for 
various percentages of mixed laminar-turbulent flow. 
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Refueling. - The program will compute a refueling segment but it must occur 
sometim '•unng cruise. Refuel is Initiated by specifying refuel at the start of 
cruise (IRFSOC = 1), or at end of cruise (IRFEOC = 1), or at a range (RNGSRF). 
During the refuel segment, the fuel taken onboard and the fuel used while cruising 
are used to find the net fuel added. Refuel continues until the gross weight of the 
aircraft with net added fuel equals the input weight after refuel (WGARF). Typical 
results for various refueling flow rates are shown in figure 34. 


CONCLUDING REMARKS 

A computer program (ASP) has been developed to allow preliminary sizing and 
performance evaluations for subsonic, supersonic, and hypersonic aircraft that use 
either JP, liquid-hydrogen, or liquid-metnane fuel. Both jet and propeller 
propulsion systems can be treated. Some of the output the program provides 
Includes: range for given gross weight and passengers, gross weight for given range 

and passengers, passengers for given range and gross weight, payload-range trade, 
passenger and fuel off-loading, effects of technology improvements, effects of in- 
flight refueling, radius to payload drop, and the basis for preparing aircraft 
sizing "thumbprints" which allow the selection of wing and engine size for best 
performance. 

The program has been used extensively for in-house studies of subsonic and 
supersonic transport aircraft concepts, some of which were developed by industry. 
For these, the results obtained by this program agreed very well with those provided 
by the companies. For several of the concepts, the program has been used to direct 
design Improvements through changes in wing and engine size. These improvements 
were later confirmed by in-depth contractual studies. 
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WEIGHT EQUATIONS 


This section contains the empirical relationships used to compute the effects of 
sizing on the weight characteristics of the aircraft. The expressions contain only 
the essential sizing parameters and, therefore, cannot be used to compute absolute 
weights. Instead, weight scaling is done by the following general form. 

W I = W I, Bl 07 f bl ) 

where Wj is the weight of the scaled item, Wj >BL , the weight of the baseline 

item, F, the result of the weight equation (s) applied to the scaled item, and f BL» 
the similar result for the baseline item. In this way, the relationships are used 
only to account for sizing changes through increments to baseline weights. At 
hypersonic speeds, aircraft require "hot-structu^e" solutions to survive severe 
local temperatures, but, as noted, the relationships do not account for this. Even 
so, the above approach still seems appropriate, at least from a preliminary 
standpoint, because the higher weights for the higher temperatures are introduced 
through the baseline aircraft. 


The weight rel ationships are listed below without proof of their validity. 
Also, since they have been obtained from several different industry contractual 
studies, a recognition of their source, in most cases, is omitted. 


Subsonic Wing Weight 


where 


K Wg N, A b (.375 + ,7X) 1R45 1 07 

F., = + .57 S + .572 S 1 * J/ 


F K S 


.-6 


K = 1.427 x 10' u (1 - 0.45 N £ )(l + /6.7'5/b) 


A 

X 

S 


number of engines on wing 
wing span 
gross weight 
ultimate load factor 

aspect ratio based on basic trapezoidal planform 
wing tip chord-root chord ratio 
wing reference area 


T8 


t/c wing thickness-to-chord ratio (weighted average) 

K s cos 2 A 

A sweep angle of the wing quarter chord line 

For a wing with multiple panels and varying t/c, the weighted average t/c is 
obtained from: 

l K t 4 (t/c) + (t/c). J 

M P r <• 

t/c = S 

b/2 

Subscripts t and r refer to the tip and root of the panel, b , the span of the 
panel, and N, the number of panels. p 


Supersonic Wing Weight 


where 


F w = .0241 


( s W 53 V 1 


T 7T 


HM 


^RM 

W RM,1 


N Z j W M ~ W RM^ 

Cos A 

“g + ") 

Vl * nacelles V " (Me * “ TR " " n) 
“f.w ! c r * *) * Vr, “s > y h f w 


An iterative solution of this cycle of equations is required until the initial 
estimate and final value of the wing weight, Fy, are about equal. Quantities 
that do not change during the iterations are: 


W MLf, = ,8803 W LG 

C R = Y D -5T^ 

* u w F,W,MAX 

If there are no fuel tanks in the wing, C R = 0. 
t hm = T a + t c 
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A I + A 0 * C I.l T 0,l * C 0,l T I,1 


= C 0,l T 0,l 
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In these equations 

S wing reference area 

b wing span 

weighted average of wing thickness- to-chord ratio 
A aspect ratio, b^/S 

A wing taper ratio 

A sweep angle of wing mid chord 

Wg gross weight 

engine weight (one) 

W TR thrust reverser weight (one) 

W^j nacelle weight (one) 

spanwise distance to nacelle-engine combination 
landing gear weight 
Wp ^ weight of fuel in wing 

Wp w maximum weight of fuel in wing 

Wg maximum width of fuselage 

The equations for both subsonic and supersonic wings were provided by Kentron 

International, Inc., Kentron Technical Center, Hampton, Virginia. 

Canard, Horizontal and Vertical Tails 

The weights of these items are directly proportional to their total wetted areas. 

Fuselage 

F f ■ (N z W s S w<c r 6 

where p Is the total wetted area of the fuselage 

Landing Gear 

f lg ■ c« l f ) -25 
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where L p * length of fuselage 

Nacelles 

The weight of the nacelles is directly proportional to their total wetted area. 

Engines 

E W 

F E ■ N E T 

where = number of engines 

T sea-level static thrust of one engine 

Ey engine weight exponent 

Thrust Keverser 

p = n T*®® 

h r« ^E 

where and T are defined above. 


F 



(.5 Ep ) 


Propeller 


where Ep propeller weight exponent. 


Gear Box 

p s n •** 

r GB V 
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Miscellaneous Systems 


The weight of this item is assumed to be unaffected by siting. 


Fuel Systems for Conventional Fuel 

Only the weight of the plumbing Is considered for these systems. For subsonic 
design Mach numbers 


Fp $ = 8.117 (N e )* 825 (T)* 38 


and for supersonic design Mach numbers 

( u • 565 

r _ l T,MAX j 

l-p_ - n 

[u+ZjM M) 

W F,MAX 


where 

Wp ^ maximum fuel load 

Wp w MAX maximum fuel load in wing 


Fuel Systems for Cryogenic Fuel 


These systems Include the weight of the plumbing, fuel tanks and insulation. 
Plumbing weight Is given by: 

F p$ = .01569 (T t )* 75 + 25.433 (l^)* 825 (T)* 38 


where 


Tj total sea-level static thrust 

The weights for the fuel tank and Its Insulation are both directly proportional to 
the total wetted area of the fuel tank. 
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Surface Controls 


If the design Mach number is subsonic 

-.5 


(Lp + '/cos a)S‘ 


SC 


7s rr 


TR 


where 


Lp = length of fuselage 
b = wing span 

A - sweep angle of quarter chord 
$ = total wing area 

A.^ = wing aspect ratio 

If design Mach number is supersonic, the weight of the surface controls is directly 
proportional to the gross weight of the aircraft. 

Auxiliary Power Unit 

The weight of this item is assumed to be unaffected by sizing. 

Instruments 

r w .354 . .5 . .68 
F IN = N E L F b 


where 


Hydraul ics 

W r .34 

: hy = ^ + Ne 


Sqj total area of surface controls 
Mq design Mach number 

Avionics 

F s u 
r A V H G 
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H metrical System 


where 


P ~ Tim 1 • ^24 /•( K -i .69 ('i i • 473 

f el = v n e J [i-p + hj i.i Av j 


W AV weight of the avionics 


Furnishings and Equipment 

The weight cf this item is directly proportional to the number of passengers. 

Air Condi tier. ing 


where 


where 


f ac ‘ ( n pas) 1-2 * n - 17 W® ( v p)’ 58 


Vp volume of passenger compartment 

N pA $ number of passengers 


f ai ■ h 


.95 


Anti -icing 


Unusable Fuel 


F up = 100 N e + .176 S 


Engine Oil 


F = u 
r E0 W E,T 


.26 


^E,T weight of all engines 


Passenger Service 

The weight of this item is directly proportional to the number of passengers. 
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Cargo Containers 


The weight of cargo containers is directly proportional to the maximum cargo 
weight. 


Crew 

Weights for the flight and cabin crew are directly proportional to the number of 
crew members. 


Passengers and Baggage 

Weights for these items are directly proportional to the number of passengers 
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DEFINITIONS OF INPUT 

The program is coded in FORTRAN Extended Version 4 (ref. 15). Provisions are 
made for two input files. One file contains propeller design data and is required 
only if propeller output data is to be computed. The data on this file must be 

input in the following format. 

Col u tin*) Type of Data 

9-10 number cf propeller blades, a right justified integer 

11-20 Mach number 

21-30 propeller advance ratio 

31-40 propeller power coefficient 

41-50 propeller thrust coefficient 

Data is listed by increasing power coefficient at constant advance ratio, then by 
increasing advance ratio at constant Mach number, then by increasing Mach number. 
During program execution, the file containing this data must be named TAPt3. 

The other file contains the input for the baseline aircraft and sizing 
options. It is always required and during execution it must be called TAPE5. The 
input segments on TAPE5 are listed below in their order of occurrence. 

Input segments for sized aircraft 

Aircraft Identification 

$PRNTC0N 

$DESGNVB 

STWINFUS 

$MISSVB 

$RESRV 

$TECHCHU 

Input segments for baseline aircraft 

Wave-Drag Geometry Deck or SSACGEm 

SAIN 

SWTIN 

SENIN 

Baseline Engine Characteristics 

Seoments prefixed by a $ sign are inputs in the NAMELIST format (ref. 15). 
Definitions of input variables follow, with default values, if they exist, enclosed 
it parenthesis behind each definition. Variables with parenthesis are arrays and 
the included numbers indicate their maximum size. 
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INPUT SEGMENTS FOR SIZED AIRCRAFT 

The first input is the aircraft identification. It is one free-field line of 
8U characters. 

r 

$PRNTC0N 

- ■} 

V-’! Inputs in this segment control program output. All inputs are integers with a 

\[ value of either 1 or 0. Unless noted, a value of 1 activates the option; 0 
suppresses it. 

\ IPPRNT = 1 prints all output 

= 0 prints only output for sized aircraft (default) 

IPAERU prints aerodynamic data, (0) 

~j IPRPRUP prints engine characteristics (U) 

3 IPRWTSM prints weight statement, (0) 

IPRMPR. r pr'nts mission profile data, (U) 

-] ISCOPE prints output showing progress of execution during interactive 

j terminal sessions, (0) 

3 

$DESGNVB 

This input segment contains aircraft design and sizing information. For 
clarity, the inputs have been separated into typical groups. 


Miscellaneous Design Input 

ISSAC = 0 for supersonic aircraft (default) 

= 1 for subsonic aircraft 

ISGEUM -- 0 if aircraft geometry is input in wave-drag format (default) 

=1 if aircraft geometry is input in $SACG£M 

UP = 1 for aircraft using JP fuel, (0) 

! IH2 = I for aircraft using hydrogen fuel, (0) 

IMTH = 1 for aircraft using methane fuel, (0) 

N0FIW = 1 for aircraft with no fuel in wing, (U) 

IBLAC * 0 for sizing of baseline aircraft (default) 

* 1 for no sizing of baseline aircraft 

RH0JP density of JP fuel, Ib/gal, (6.75) 
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RH0H2 


density of liquid hydrogen fuel, lb/gal, (.592) 


OtiLW 

NPPUGL 

STSAB 

LEVELS 


FHDES 

IPASDES 

NENDES 

NPOES 

RMAXOES 

ZLBDES 

NCARDES 

TUWDES 

TPEDES 


width of door-galley lavatory, ft, (7.3) 

number of passengers per door-galley-lavatory, integer, (75) 

number of seats abreast; 0 provides the baseline value STSABR, (0) 

number of passenger levels, 0, 1, or 2; 0 provides the baseline 
value, (0) 

wall thickness of passenger cabin, inches, (4) 

wall thickness of fuel tanks for liquid hydrogen or methane fuel, 
inches (9) 

Specific Design Inputs 
design fuselage fineness ratio 
number of passengers, integer 
number of engines, integer 
number of engine nacelles, integer 
design radius of fuselage, ft 
design length of fuselage, ft 
design cargo load, lb 
design th.'ust-to-weight ratio 
design single engine thrust, lb 
design wing loading, psf 



RHOMTH 

density 

of liquid methane fuel, lb/gal, (3,54) 

i5 

i 

A 


RHOCAR 

density 

of cargo, 1 o/ft (10. ) 

| 

t 

j 


EIPEN 

reduction in thrust to account for interference effects between 
propeller and wing, percent, (0) 

5 


The 

next three variables apply to aircraft using hydrogen or methane fuel. 

j 

\ 

* 

KBFN 

- 1 

for constant fuselage fineness ratio, (0) 

i 

j 


KBL 

= 1 

for constant fuselage length, (1) 

j 

t 

0 

KBR 

= 1 

for constant fuselage radius, (0) 
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Passenger Cabin Design Input 





WOSDES 


REFAUES 





design wing reference area, ft^ 

WGOES design gross weight, lb 

RNGUES design range, n.mi. 

ZMCROES design cruise Mach number 

The next six variables apply to subsonic aircraft. 

ASPROES design wing aspect ratio 

EDES design endurance, h 

POWOES design shaft horsepower-to-wei ght ratio, hp/lb 

SHPDES design shaft horsepower, hp 

SQCUES design sweep of wint> quarter chord line, deg 

WTOCUES design thickness ratio of wing, percent 

If RNGUES is not input, the program solves for the range at input gross weight and 
passenger load. If RNGUES is input, the program so-ves for the gross weight that 
provides this range at input passenger load. 

To find the passenger load for a given range and gross weight, use the following 
integer variable. 

IGRAWFP = 1 (default is 0) 

along with inputs for WGDES and RNGUES. 

Each of the above variables, which defaults to its baseline value, is used to 
hold that design parameter constant through one or more sizing passes during program 
execution. If design parameters are to be changed in each sizing pass during 
execution, the following design variable arrays must be used. All arrays are one- 
dimensional. 

Ue sign Variable Arrays. - 

TOWTB (7 ) an array containing design thrust-to-wei ght ratios 

TPETB ( 7 ) an array containing design thrusts for a single engine, lb 

NTOW the number of values in either the TUWTB or TPETB arrays, (0) 

W0STB(7) an array containing design wing loadings, psf 

REFATb ( 7 ) an array containing design wing reference areas, ft^ 

NWOS the number of points in either the WOSTB or REFATB arrays, (0) 

ARTrt(lU) an array containing design wing aspect ratios 
30 




NASPR 


the number of points in ARTB, (0) 

POWTtJ ( 7 ) an array containing shaft horsepowe --to-weight ratios 
NPOW the number of values in POWTB, (0) 

Restrictions on Use of Variables .- Not all of the design variables are 
independent parameters, the following restrictions, therefore, are placed on the 
use of several of these variables. 

1. Since POWUES, TOWDES, TPEOES, POWTB, TOWTB, and TPETB all control engine 
sizing, use only one of these variables. 

2. Similarly, WOSDES, REFAUES, WOSTB, and REFATB all affect wing sizing, 
therefore, use only one. 

3. The ARTS array must not be used with any other design variable array. 

4. The integer variable I6RAWFP may not be used with any design variable array. 


$TW1NFUS 

This input segment contains input for the geometry of twin-fuselage aircraft. 
The only component of this aircraft that can be sized is the propulsion system; 
therefore, only the design variables WGOES, TOWDES, TPEOES, and RNGOES are 
allowed. If changes in any other design variable, like wing size, passenger load, 
etc., are required, a new baseline aircraft (with its aerodynamics weights and 
propulsion) must be developed and input. 


The only geometry input necessary for twin-fuselage aircraft are the reference 
lengths and wetted areas of the components. These are required for skin-friction 
calculations during the mission analysis. To prevent program aborts, however, a 
complete geometry input segment in either the wave-drag format or $SACGEM is also 
required but any available input segment will suffice because it will iot affect 
results. 


ITWINF 

FUKUN 

FUSWTIN 

TWFREFA 

IWSEGIN 

WRLIN(2U) 


= U for conventional aircraft, (default) 

= 1 for twin-fuselage aircraft 

fuselage reference length, ft 

fuselage wetted area, total of both fuselage, ft 

wing reference area, ft 1 * 

number of wing segments 

a one-dimensional array containing the reference length of each wing 
segment, ft 


WSWETIN(20) a one -dimensional array containing the wetted area of each wing 
segment, ft^ 
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CRL IN canard reference length, ft 

CSWETIN canard wetted area, ft^ 

HRLIN horizontal tail reference length, ft 

HSWETIN horizontal tail wetted area, ft^ 

IFINS number of vertical fins 

FRL IN (6 ) a one-dimensional array continuing the reference lengths of the vertical 
fins, ft 

FSWETIN(o) a one-dimensional array containing the wetted areas of the vertical 
fins, ft^ 

PRLIN the reference length of one engine pod, ft 

PSWETIN the wetted area of one engine pod, ft^ 


SMI SSVS 

This segment provides the inputs requir'd for the flight profile analysis, 
inputs, again, have been separated into typical groups. 


The 


Overall h.ssion Input 

IENOUR = 1 mission profile for long endurance aircraft, integer, (0) 

I RADIUS = 1 provides range to payload drop, integer (0) 

DEW operating empty weight, lb, (if greater than U, this value will replace 

the calculated OEW) 

PAYLOAD total payload weight, lb, (if greater than 0., this value will replace 
the calculated payload) 

POLF payload factor, percent (100% gives full payload), (100.) 

NPASOL number of passengers in off-loaded aircraft, integer (-1, the default, 
gives no off-loading) 

WCAROL cargo weight in off-loaded aircraft, lb 
(-10, the default, gives no off-loading) 

ICALPRC = 1 provides data for the range-payload curve, integer, (0) 

OLRTB(IO) a one-dimensional array containing the ranges for fuel off-loading, 
n.mi. 


NOLR 


the number of values in OLRTB, integer, (0) 
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RRTB(IO) 

NRR 

IRFSOC 

IRFEOC 

RNGSRF 


RFRIGPM 

WGARF 

UELTCG 

EM 

TDFLM 

NOTO 

FWOTO 

ALFDOT 

DELTCTO 

OTGRUP 

I1EOUT 

RFF 

TXTME 

GRALFA 

TOALFA 

TOOELF 

TOFLTO 


a one -dimensional array containing reduced ranges with no fuel 
off-loading, n.mi. 

the number of values in RRT6, integer, (0) 

= 1 start refueling at start of cruise, integer, (0) 

= 1 start refueling at end of cruise, integer, (0) 

range to start of refuel, n.mi. 

(refueling must occur sometime during cruise. If RNGSRF is set too low, 
refuel will begin at the start of cruise; if set too high, refuel will 
begin at the end of cruise. To activate refueling, the value of RNGSRF 
must be greater than 10.), (0) 

refuel flow rate, gal /min, (600.) 

the aircraft gross weight after refueling, lb 

(-10, the default value, gives the design gross weight) 

standard day temperature increment during main mission, deg., C (0.) 

surface emissivity, (.8) 

thrust deflection during main mission, deg., (0) 


Takeoff Input 

= 1 no takeoff calculations, (0) 

fuel weight used during takeoff, lb 

(may also be input as a fraction of gross weight), (0.) 

aircraft rotational speed, deg/sec, (3.) 

standard day temperature increment during takeoff, deg, C, (0. ) 
time for landing gear retraction, sec, (10.) 

* 1 one engine out during takeoff, integer, (0) 
rolling friction factor, (.02) 
taxi time, min, (10.) 

angle of attack before start of rotation, deg (-4) 
angle of attack after rotation, deg, (5.5) 
flap deflection during takeoff, deg, (30.) 
thrust deflection during takeoff, deg, (0.) 
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TFACT 


factor for increasing or decreasing net thrust during takeoff, (1.) 
WFFACT corresponding factor for fuel flow, (1.) 

HOBS obstacle height, ft, (35.) 

HNOGE altitude for disappearance of aerodynamic in-ground effect, ft, 

(wing span) 

TOCLGR takeoff climb gradient, percent, (6.8) 

Climb Input 

HA$NT(50) a one-dimensional array containing the altitudes in the climb profile 

ZMASNT(50) a one-dimensional array containing the corresponding Mach numbers in the 

climb profile 

ICLPR the number of points in the climb profile 

ZMSTTH Mach number at which engine throttling may begin to conserve climb fuel, 

(. 6 ) 

ZMMPCL Mach number above which maximum engine thrust is used, (1.) 

Cruise Input 

NCRP the number of points at which cruise calculations are made, (4) 

ZMCR cruise Mach number 

ROCMIN the minimum rate of climb during cruise, ft/min 

CRALT1 altitude for start of cruise, ft ( - 1 U ) 

CRALT2 altitude for second leg of two-step cruise (-10) 

lbFCR = 1 gives cruise at best available lireguet factor, integer, (0) 

ICACR = 1 gives cruise at constant altitude, integer, (0) 

ICLCR = 1 gives cruise at constant lift coefff icient, integer, (0) 

ALFCR angle of attack for cruise, deg 

(Note: this input may be used only if CLAR or ALPHAT are input in 

NAMELIST segment $AIN) 

CRTOTM cruise throttling factor (0.) 

(if CRTOTM = 1., full power is applied; if defaulted, CRTOTM is not 
used) 
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The following input is for subsonic cruise legs on a supersonic mission. Range 
or time for cruise may be input, but not both at the same time. Cruise on both legs 
may occur. 

For subsonic cruise on the outbound leg: 

ZMSBCOB cruise Mach number 

HSBCOB cruise altitude, ft 

SSRNGOB cruise range, n.mi. 

TSBCOB cruise time, min 

For subsonic cruise on the inbound leg: 

ZMSBCIB cruise Mach number 

HSBCIB cruise altitude, ft 

SSRNGIB cruise range, n.mi. 

TSBCIB cruise time, min 

Input for Long-Endurance Cruise 
CLMAX maximum lift coefficient, (2.2) 

NENCR number of engines operating during cruise (-1) 

DCFACT increment in drag coefficient for engines not operating, (0.) 

VKCR1 velocity at start of cruise, k 

ICONVCR = 1 for constant cruise velocity, (U) 

Descent Input 

ICALUS = 1 descent calculations are made 

= 0 descent calculations are not made 

NUSP the number of points at which descent calculations are made, (4) 

USFF the fraction of fuel flow at maximum power that is used during descent, 

(.067) 

WFURCF an initial estimate of the fuel used during descent normalized by 

aircraft gross weight, (.005) 
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Approach Input 

APALFA angle of attack during approach, deg, (-40) 

APDELF flap deflection during approach, deg, (30.) 

IFXAFLP = 1 uses input flap deflection 

= 0 finds flap deflection for minimum drag, (default) 

(used only if VAPIN is input) 

VAPIN approach speed, k, (-10) 

(if defaulted, approach speed is calculated internally) 

Note: VAPIN is used only if aerodynamic data for more than one flap 

deflection is input in $AIN. 

$RESRV 

The input for the reserve fuel legs are placed in this segment. 

RFWT reserve fuel weight, in either pounds or as fraction of gross weight. 

(use this input only if reserve fuel calculations are to the bypassed) 

IRCRO = 0 calculates climb, cruise, and descent segments of reserve 

fuel legs 

= 1 calculates only cruise segments of reserve fuel legs, (!' 


TPFATME 


ICALMA 


RMAFA 


Trip Fuel Allowance 

time at fuel flow at start of cruise, h or min 
(use for subsonic aircraft), (0.) 

percent of trip fuel 

(use for supersonic aircraft), (0. ) 


Missed Approach 

= 0 uses RMAFA to calculate fuel requirement 
= 1 calculates missed approach to find fuel requirement, (1) 

time at takeoff fuel flow, min, (2.0) 


Alternate Airport 

AAU distance to alternate airport n.mi., (0.) 

For this leg, the program normally iterates through Mach number, altitude, and 
engine thottling during climb to find the conditions for least fuel. If any of the 
following variables are input, iteration on that condition will not be performed. 

HALT altitude for cruise to alternate airport 
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TUTMALT factor for engine throttling during climb to cruise, value between U 
and 1. 


ZMALT 


Mach number for cruise to alternate airport 


time for hold, .nin 
altitude for hold, ft 


ZMHLD 


Mach number for hold 


Again, optimum conditions for least fuel are normally determined but inputs in HHLD 
and/or ZMHLD will eliminate the iteration on that condition. 

$TECHCHG 

This segment contains the inputs for finding the effects of changes in the 
technology of the various aircraft systems. 

TIML the lowest Mach number at which changes apply, (-1) 

TIMH the highest Mach number at which changes apply, (10.) 

TDCUT (15) a one-dimensional array of percentage changes in total drag coefficient, 
a function of Mach number; a positive value gives drag reductions 

TMDCUT (15) the Mach numbers corresponding to the values in TOCUT 


NDCDT 


the number of values in the TUCDT array, (U) 


The following innut is all in percentages. A positive value gives a technology 
improvement. Default values for all inputs are 0. 

DSFC engine specific fuel consumption 


DCDW wave drag 

DCURUF roughness drag 

DCUAC air-conditioning drag 

DCDBL bleed drag 

DCDO zero-lift drag 

DCDL drag due to lift 

DCDT total drag 


T'H m 


- r 



Friction Drag Changes 


DCDFWNG 

DCUFCS 

DCDFBUY 

UCDFPDS 

DCDF 

DWWING 

DWCHV 

DWFUS 

DWLG 

UWNAC 

DWST 

DWENG 

DWTR 

DWTPR 

UWTGB 

L)WM I S 

DWPLM 

DWFT 

DM I NS 

UWPPS 

DWSC 

DWAP 

DM INST 

DWHYD 


wing 

canards, horizontal tails, and vertical fins 

fuselage 

nacelles 

total friction drag 


Weight Changes 

wing 

canard, horizontal tail, vertical fin 

fuselage 

landing gear 

nacel les 

total structure weight 
engi nes 

thrust reversers 
propel lers 
gear box 

miscellaneous propulsion system weights 
plumbing system 

fuel tank for hydrogen or methane fuel 

fuel tank insulation for hydrogen or methane fuel 

total propulsion system weight 

surface controls 

auxi liary power 

i nstruinents 

hydraulics 

electrical 


DUELE 


A 




DWAVON 

avionics 

DWFEO 

furnishings and equipment 

0WAC0N0 

air conditioning 

OWANTICE 

anti -icing 

OWSEQ 

total systems and equipment weight 

OWEMP 

empty weight 

DWOEW 

operating empty weight 


Laminar Flow 

The following inputs express the percentage of laminar flow on the surface. A 
value of 0. gives all turbulent flow; a value of 10U. gives all laminar flow. All 
default values are 0. 

PCLAMW wing 

PCLAMF fuselage 

PCLAMC canard 

PCLAMH Horizontal tail 

PCLAMBV body vertical fins 

PCLAMWV wing vertical fins 

PC LAMP nacelles 


INPUT FOR BASELINE AIRCRAFT GEOMETRY 

The geometry for the baseline aircraft can be input by one of two methods. The 
first method which applies only to supersonic aircraft uses the zero-lift wave-drag 
program format described in reference 7. This input is very detailed and unless it 
is already available from the analysis of wave drag, the effort required to assemble 
the information is not justified for this program. The alternate method available 
is the use of $SACG£M. 


SSACGEM 

Wing 

REFA wing reference area, ft^ 

WXO distance from fuselage apex to wing apex at wing centerline, ft 
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WRC 

root chord at wing centerline, ft 


WTC 

tip chord, ft 


WSPAN 

wing span, ft 


WDXT 

longitudinal distance from wing apex to leading edge of tip chord, ft 
(positive rearward) 

YWCREX 

spanwise distance from wing centerline t.o exposed root chord, ft 


NWBP 

number of breakpoints in wing leading and trailing edges, an 
Wing apex and tip are not included. If leading edge and trail 
break at the same spanwise station, count as one breakpoint. 

i nteger. 
ing edge 

WDX8P ( 12 ) 

a one-dimensional array of the longitudinal distances from wing 
leading edge at breakpoints, ft (positive rearward) 

apex to 

WYBP(12) 

a one-di Tiensi onal array of the spanwise distances to the breakpoints ft 

W8PC (12 ) 

a one-dimensional array of wing chords at the breakpoints, ft 


WTOC 

average maximum wing thickness-chord ratio in percent 



Fuselage 


FLGTH 

length, ft 


FDPTHMX 

maximum depth, ft 


FWDTHMX 

maximum width, ft 


FRAD 

maximum radius, ft 

(Note: If FRAD is input, FDPTHMX and FWDTHMX inputs are not required) 

FWETA 

total fuselage surface area, ft/ 


FVULTUT 

total fuselage volume, ft 



Horizontal Surfaces 


The dimensions tor horizontal surface*, are input in the toll oi wng one- 
dimensional arrays. 

NH 

the number of horizontal surfaces, (2 maximum) 


HDXT (2 ) 

the longitudinal distance from the apex of the surface to its 
edge at the tip, it 

leading 

HRC (2 ) 

root chord, ft 


HTC (2 ) 

tip chord, ft 
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HSPAN (2 } 
HTOC (2 ) 
HXO (2 ) 


span, ft 

average thickness- to -chord ratio, percent 

longitudinal distance from fuselage apex to surface apex, ft 


Vertical Surfaces 

The dimensions for the vertical surfaces are input in the following one- 
dimensional arrays. A spanwise mounted surface is assumed to be located 
symmetrically about the fuselage centerline and is counted as one surface. 

NV the number of vertical surfaces, (4 maximum) 

VUXT(4) the longitudinal distance from the apex of the surface to its leading 

edge at the tip, ft 

VRC (4 ) root chord, ft 

VTC(4) tip chord, ft 

VSPAN (4 ) span, ft 

VT0C(4) average thickness-to-chord ratio, percent 

VXO (4 ) longitudinal distance from fuselage apex to surface apex, ft 

VYO ( 4 ) spanwise distance from fuselage centerline to surface, ft 

Nacel les 

X0P(9) a one-dimensional array containing the longitudinal distances from the 
fuselage apex to the nacelle, ft 

YOP (9 ) a one-dimensional array containing the spanwise distances from the 

fuselage centerline to the nacelle centerline, ft 

NPI the number of nacelles in the above arrays 

Spanwise mounted nacelles are counted in the same manner as the vertical 
surfaces. 

$A1N 

The baseline aerodynamic input is located in this segment. Trimmed angle of 
attack, lift, and drag data are preferred. 
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TALPTUG(I5,4) 

T | 'LT0G(15,4) 
TCDTOG ( 15,4 ) 
NTOG 


TALPHTU(15,4) 
TCLT0(15,4) 
TCDTO ( 15,4) 
NTOP 

NFU 

TFSET (4) 

CDGRT (15) 

CLGKT (15) 

NLGU 

DCOLG 


CLT (15,15) 

ALPHAT (15,15) 
CDPT (15,15) 
NAI 


Takeoff Aerodynamic Input With Ground Effect 

a two- ii mens ional array of angle of attack for various flap 
deflections, deg 

a similar array for the lift coefficients, 
a similar array for the drag coefficients, Cq 

the number of values in each of the above arrays at each flap 
deflection, (15 maximum) 

(If NTOG = 0, an internal routine is used to calculate values in the 
above arrays when the following takeoff input is provided.) 


Takeoff Aerodynamic Input Without Ground Effect 
angle-of -attack array 
lift-coefficient array 
drag-coefficient array 

the number of values in the above arrays at each flap deflection, 

(15 maximum), (U) 

the number of flap deflections, (4 maximum) (0) 

a one-dimensional array containing the flap deflections, deg 

a one -dimensional array of increments in drag coefficient of the 
landing gear as a function of lift coefficient 

the corresponding array of lift coefficients 

the number points in the above arrays, (0) 

increment in drag coefficient of landing gear 
(independent of lift coefficient, use if NLGU = U 


Main Mission Aerodynamic Input 

a two-dimensional array of lift coefficients for various Mach 
numbers 

the nrresponding array for angles of attack, deg 

the corresponding array for the drag coefficients 

the number of values in each of the above arrays at each Mach number 
(the value of the first dimension; NAI points must bo input at each 
Mach number) 
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IPBOL = 0 aerodynamic input in CLT, ALPHAT, and COPT will be used for 

aerodynamic data (the default value) 

= 1 the parabolic polar inputs (CDOR, CLAR, OLFR, and CLOR 
defined below will be used for aerodynamic data 

CD0R(15) a one-dimensional array containing the minimum drag coefficients at 

various Mach numbers 


CLAR ( 15 ) 
DLFR(15) 
CLOR (15) 
NAJ 

MAER0T(15) 


a similar array for the lift-curve slope, per deg 

the array for drag-due-to-1 ift factors 

the array for lift coefficients at minimum drag coefficient 

the number of Mach numbers for the lift and drag data (15 maximum) 

a one-dimensional array containing the Mach numbers for the lift and 
drag data 


THAR0(15) a one-dimensional array containing the altitudes at each Mach number 

at which the skin friction contributions to drag were calculated 


DELC0(15) a one-dimensional array containing increments in zero-lift drag 

coefficient for the Mach numbers in MAEROT 


COW(IS) 

C0RUF(15) 

C0AC(15) 

C0BL(15) 

Note: 
or CDOR. 


a one-dimensional array containing the wave-drag coefficients at 
each Mach number in MAEROT 

a similar array of roughness drag coefficients 

a similar array of air-conditioning drag coefficients 

a similar array of engine-bleed drag coefficients 

Input for COAC and CDBL is required only if values are included in CDPT 


$WT1N 

This segment contains the weights and all other input for the baseline aircraft 
but not the engine data. Except where noted default values are 0 and input weights 
are in pounds. 


Structural Weight Input 


WWING wing 

WHT horizontal tail 
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WVT 

WCAN 

WFUS 

WL(i 


WFT 
WINS 
WFT INS 

TSA 

WSC 

WAP 

WINST 

WHYD 

UELE 

WAVON 

WFEQ 

VJACOND 

WANT ICE 

WFCP 

IFCR 

WCCR 

ICCR 

WUFlIEL 





vertical tails (sum of all) 

canard 

fuselage 

landing gear 


Weight Input for Cryogenic Tanks 
(Hydrogen on Methane Fuel ) 

fuel tank, total 
insulation, total 

combined weight of fuel tank and insulation (required only if WFT and 
WINS cannot be separated) 

total surface area of fuel tank, ft 

System and Equipment Weight Input 
surface controls 
auxiliary power 
i nstruments 
hydraul ics 
electrical 
avionics 

furnishings and equipment 
air conditioning 
anti -icing 


Operating Weight Inputs 
weight of flight crew 
number of flight crew members, integer 
weight of cabin crew 
number of cabin crew members, integer 
unusable fuel 
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WE NO 
WPSER 
WCCONT 
WOPIN 


WPAS 

IPASR 

WPB 

WCAR 


WGKEF 

ULF 

ZMDESBL 

STSABR 

LVLR 

XPASC 

VCARMX 

BFCILBS 

WFCILBS 

BFCIGAL 

WFCIGAL 


engine oil 
passenger service 
cargo containers 

sum of all operating weights (use only if all individual weights cannot 
be supplied) 

Payload Weight Input 

weight of passengers 
number of passengers, integer 
passenger baggage 
weight of cargo 


Miscellaneous Input 
gross weight of baseline aircraft 
design load factor, (3.5) 
design Mach number for baseline aircraft 

number of seats abreast 1 

number of passenger levels, (2 maximum) i 

length of passenger compartment, ft 
maximum volume for cargo, ft 

maximum fuel capacity in fuselage, lb „ 

maximum fuel capacity in wing, lb 
fuel capacity in fuselage, gal 
fuel capacity in wing, gal 

(Note: input for I'FCIGAL and WFCIGAL is not required if BFCILBS and 

WFCILBS are used) 

i 

t 

I 






^ r-t 


$ENIN 

Data for the engine and nacelle are input in this segment. 


NENR 

the number of engines, integer 


NEW 

the number of engines mounted on the wing, integer 

NEF 

the number of engines mounted on the fuselage. 

integer 

NPODSR 

the number of engine nacelles, integer 


PDAVG 

the average diameter of the nacelle, ft 


PLGTH 

the length of the nacelle, ft 


ESEXP 

nacelle length sizing exponent, (.438) 


EWEXP 

engine weight sizing exponent, (1.085) 


SLSTHR 

sea level static thrust for the engine having 

the following weight, lb 

WENS 

weight of one engine, lb 


WTR 

weight of one thrust reverser, lb 


WNAC 

weight of one engine nacelle, lb 


WMIS 

total miscellaneous propulsion system weight. 

lb 

WPLM 

total propulsion plumbing system weight, lb 



The following three variables are provided for engine bleed and air-condition 
drag. These variables should be used only if the drag items are not included in 
gross thrust values or by entries in CUAC and COBL in $AIN. 


D(X|1NT(15) a one-dimensional array containing the sum of engine bleed and air- 
condition drags divided by dynamic pressure for a range of Mach numbers, 


DUQMT(15) the corresponding arrav of Mach numbers 
NBDP the number of Mach numbers in the above arrays 

Note: These variables cannot be used for engines with two-dimensional inlets. 


input for Propeller-Engine Combinations 


Propeller-driven aircraft require no additional input if: (1) propeller output 
data is not required, (2) the weight of the propeller and gear box are included in 
the engine weight, WEN6, and (3) the design inputs TOWUES, TPEUES, TOWTB, or TPETB, 
rather than PUWDES, SHPDES, or POWTB are used to size the propulsion system. This 
analysis treats the aircraft as a jet -driven aircraft. 
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The propeller wilt be treated in more detail, however, if TAPE3, described 
previously, is provided and the following data is input. 

AUVRUES design advanced ratio 

CPUES design power coefficient 

CTOES design thrust coefficient 

EPRDES propeller efficiency at design conditioning 

EOVDES overall efficiency of the propeller-engine combination at design 
conditions 

P0D2DEi> shaft horsepower-propeller diameter squared ratio at design conditions, 
hp/ft* 

PUPS propeller tips speea at design ronditons, fps 

DIAMP8L propeller diameter, ft (-10.) 

SHPREF shaft horsepower of the engine at sea level static conditions and having 

its weight entered in WENG 

RPMEBL engine RPM 

RPMPBL propeller RPM 

PSEXP propeller weight sizing exponent (2.488) 

WTG8 weight of one uearbox, lb (-10.) 

WTPR weight of one propeller, lb (-10.) 

The propeller weight, WTPR, is used to test if a jet analysis the (default case) 
or a propeller analysis is to be made. When WTPR is input, all of the above 
propeller inputs and the P/U^ values, described below, are required. 


Input for Engine with Two-Dimensional Inlets 

1NL2U = 1 engine nacelles are two-dimensional 

= 0 engine nacelles are circular, (default) 

PWIU2D width of the nacelle, ft 

PHGT2D height of the nacelle, ft 

PLGTH length of the nacelle, ft 

wetted area of one nacelle, ft* 


PSWET2D 



BASELINE ENGINE CHARACTERISTICS 


This final Input segment contains four groups of data for the Installed 
engine. These Include the identification of the baseline engine, full-power data as 
a function of altitude and Mach number, part-power data, and for engines with two- 
dimensional Inlets, full-power data as a function of Mach number, angle of attack, 
and altitude. The engine identification, the first input of this segment, is one 
line of free-field input limited to 80 characters. 

The next input provides the full -power characteristics of the engine as a 
function of Mach number and altitude. The data format is as follows: 


Col unm 

Type of Data 

1-5 

Mach number 

6-15 

altitude, ft 

21-30 

gross thrust, lb 

31-40 

ram drag, lb 

41-50 

tuel flow, lb/h 

51-60 

P/D 2 , hp/ft 2 


Where P/D 2 is the shaft horsepower-propeller diameter squared ratio. P/D 2 data 
is required only if the propeller data described above is input. 

The data is arranged by increasing altitude at constant Mach number, then by 
Increasing Mach number. Up to 15 separate values of altitude and 15 separate values 
of Mach number are allowed. On the last line of this data segment, the characters 
9. located in columns 2 and 3 are required. 

The next group of data pertains only to engines with two-dimensional inlets and 
Is required only if the effects of angle of attack are to be included in the full- 
power data. To Indicate the presence of this data group, which is restricted to 
supersonic and higher Mach numbers, the first line of Input must contain the three 
characters -9. In columns 2, 3, and 4. On succeeding lines thereafter, the data is 


Input In the following format. 
Col umn 

Type ot ,'atd 


1-5 

Mach number 


6-15 

altitude, ft 

♦ 

16-20 

angle of attack, deg 


21-30 

gross thrust, lb 

* 

31-40 

ram drag, lb 


41-50 

fuel flow, lb/h 



The data Is arranged by Increasing altitude at constant angle of attack and 
constant Mach number, then by increasing angle of attack at constant Mach number, 
then by Increasing Mach number. Input Is restricted to 15 values of Mach number, 3 
values of angle of attack at each Mach number, and 3 values of altitude at each 
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angle of attack— Mach number combination. Sufficient data must be input to provide 
the same number of angles of attack at each Mach number and the same number of 
altitudes at each angle of attack— Mach number combination. If this requirement is 
not met, the program will abort. To end this data group, the two characters 9. must 
be located in columns 2 and 3 on the last line of input. 

The last input provides the part-power characteristics of the engine. The 
format for this data is: 


Column 

Type of Data 

1-5 

Mach number 

21-30 

gross thrust, lb 

31-40 

ram drag, lb 

41-50 

fuel flow, lb/h 


Data is arranged by increasing thrust at constant Mach number, then by increasing 
Mach number. The last input of this group must contain the two characters 9. 
located in columns 2 and 3. 


SAMPLE INPUT LISTING 

An input listing that illustrates the content of the previous section is given 
in Table Bl. This input is for the design point of the supersonic transport (JP 
fuel) shown in figure 22 and the resulting output is discussed in Appendix C. 

While the sizing "thumbprint" (fig. 22) provides fairly accurate values of W/S 
and T/W for the design point, the gross weight value is very approximate. To more 
accurately determine this value, the computer must be directed to find the gross 
weight that provides the design range at the design point conditions. This is done 
in input segment $UESGNVB by inputs for the required range (KNGDES = 4000.), wing 
loading (W0S0ES = 82.), and thrust -wei ght ratio (T0W0ES = .32). These and the other 
inputs for the sized aircraft are shown in Table Bl. 


ENDIX C 


OUTPUT 

The output from the program Is contained on three mass storage files. Unless 
clanged by the user at execution time, these files are labeled TAPE6, TAPE11, and 
TAPE14. TAPE 6 contains a listing of all the results provided by the program, 
TAPE11 contains selected data In a format suitable for Input to a plotting program, 
and TAPE14 contains the basis for preparing sizing "thumbprints," 


TAPE6 

On this file, several types of output are available. These include the 
preliminary output, sized aircraft output, a short listing of selected output 
quantities, weight statement, and a mission profile output. 

Prel 1 mi nary Output. - This printout, which is optional and obtained with 
IPPRMTH in $PRl1YC0N, contains a listing of the input values for all input 
segments. Because of its length, an example of this output is not given here. 

Si/en Aircraft output .- This printout, which cannot be suppressed, contains the 
results for each sizing variation requested, A listing of the output that results 
from the input for the design-point aircraft given in Table B1 is given in 
Table Cl. A short summary of this listing is also output and this is shown in 

Table C2. A listing of the weight statement, if selected by IPRWTSM=1, is given iri 
Table C3. With IPRMPRF-1, mission profile results are output and these are listed 
in Table C4. 


TAPE11 

This file, produced for every computer cycle, contains the same parameters as 
the short output list (Table C2) located on TAPE6. The file is formatted to allow 
the sizing and performance results to be plotted during interactive plotting 
sessions where the abscissa and ordinate of the plots are identi f ied by the indices 
of the chosen variables. The Indices are the numbers of the variables on the short 
output list in Table C2. An illustration of the listing this file contains is shown 
in Table C5 and the different sections of data correspond to different computer 
cycles in which the sizing parameters are varied. The first section contains the 
same results as the short oucput list (Table C2) for the design point aircraft; the 
remaining sections contain the results for a systematic variation in thrust-to- 
weight ratio at constant wing loading. The Integer at the beginning of each section 
has a constant value for all sections of data that belong to the same curve (in this 
case, a constant value of wing loading). 


TAPE14 

This file contains the basis for preparing sizing "thumbprints" like that shown 
in figure 22. Unlike TAPE6 and TAPE11, this file Is produced only when a matrix of 
wing and engine size variations are to be computed. To produce the "thumbprint" for 
the supersonic transport, (fig. 22), inputs were required for the design range 
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(RNGOES - 4000. ), an array of values for wing loading (WOSTB = 110., 100., 90., 80., 
70., 60., 50.), an array of tnrust-to-weight ratios (T0WTB = .25, .30, .35, .40, 
.45, .50, .55), along with NW0S - 7 and NT0W = 7. The values that resulted from one 
computer run are shown in Table C6. The listing contains the abscissa values of 
wing loading, the ordinate values of thrust-to-weight ratio, and the values of the 
parameters to be contoured that are functions of wing and engine size. The rows 
provide the change in data with W/S at constant T/W, whereas the columns give the 
data change with T/W at constant W/S. The large negative values (-1000000) are 
default values used where no solutions were obtained. 

Table C6 identifies all the parameters that are available for contouring. To 
obtain the coordinates of the actual contours, the data on this file can be either 
cross-plotted by hand or TAPE14 can be used to prepare an input file for a 
contouring program. 



APPENDIX D 
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JOB CONTROL CAROS 

The program executes on the Network Operating System (NOS) Version 1.4 currently 
in use at the Langley Research Center. The control cards required for execution 
are: 

GET, ASP/UN = 273347N 
ASP (F5, F6, Fll, F14) 

where 

F5 is the file containing the baseline aircraft input and sizing 
options (TAPES) 

F6 is the output file (TAPES) 


Fll is the file containing plot data (TAPE 11) 

F14 is the file containing the data for contour plots (TAPE14) 


Note: Although the file containing the propeller design data (TAPE3) does not i 

appear in the execution card ASP, it must be a local file du r ing execution. : 

I 

The storage requirement for program execution is about 21U,U00 (octal) words. j 
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(d) Supersonic transport. Liquid hydrogen fuel. 
Figure 1. - Continued. 
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(c) Alternate airport, 260 nm (d) Hold, 30 min 

Figure 11.- Effect of sizing on reserve fuel. Supersonic transport 
JP fuel, W r = 718,000 lb 
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Figure 12 Effect of cruise options on altitude and range. 

Supersonic transport, JP fuel, T/W = .367, 

W G = 718,000 lb 
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Figure 16 Excess thrust during c ! fmb. Supersonic 
transport JP fuel, W q- 71 8,000 !b 
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Figure 27,- Effect of sizing on mission radius. Supersonic 
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